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Csl ■ Abstract 

Constraints are found on the spatial variation of finite-time Lyapunov ex- 
ponents of two- and three-dimensional systems of ordinary differential equa- 
tions. In a chaotic system, finite-time Lyapunov exponents describe the aver- 
age rate of separation, along characteristic directions, of neighboring trajecto- 
ries. The solution of the equations is a coordinate transformation that takes 
G I initial conditions (the Lagrangian coordinates) to the state of the system 

at a later time (the Eulerian coordinates). This coordinate transformation 
naturally defines a metric tensor, from which the Lyapunov exponents and 
^Sj ■ characteristic directions are obtained. By requiring that the Riemann curva- 

K>" . ture tensor vanish for the metric tensor (a basic result of differential geometry 

in a flat space), differential constraints relating the finite-time Lyapunov ex- 
^P I ponents to the characteristic directions are derived. These constraints are 

realized with exponential accuracy in time. A consequence of the relations is 
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<^ • that the finite-time Lyapunov exponents are locally small in regions where the 

^~^ ■ curvature of the stable manifold is large, which has implications for the effi- 

ciency of chaotic mixing in the advection-diffusion equation. The constraints 
^ I also modify previous estimates of the asymptotic growth rates of quantities 

^ ' in the dynamo problem, such as the magnitude of the induced current. 
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In a dynamical system, the average rate of stretching of infinitesimal elements 
of phase-space is described by the Lyapunov exponents. For a fluid, this stretch- 
ing is crucial to the enhanced diff"usion observed in chaotic flows. Though the 
exponents converge to a time-asymptotic limit that is independent of the initial 
condition, in most situations of interest this convergence is very slow: there exist 
regions where, for finite times, the exponents remain anomalously small. Infor- 
mation on the spatial and temporal behavior of finite-time Lyapunov exponents 
is thus useful in identifying such regions of low stretching rates. Our approach 
is to regard the time evolution of the solution to a set of ordinary diff"erential 
equations as a coordinate transformation. For smooth dynamical systems, we 



can study this transformation with the tools of differential geometry. Specifi- 
cally, we consider the Riemann curvature tensor, a quantity that must vanish 
in every coordinate system when using the usual Euclidean distance. By exam- 
ining the details of how the different terms of the Riemann tensor balance each 
other, we obtain constraints on the finite-time Lyapunov exponents, clarifying 
and generalizing previous work.Ercl These constraints take the form of differential 
relations between the exponents and the characteristic directions of stretching 
in the flow. The constraints have implications for chaotic mixing, where they 
allow identiflcation of locally small finite-time Lyapunov exponents, and for the 
kinematic dynamo problem, where they imply a slower growth rate of the power 
needed to sustain a dynamo. 

I. INTRODUCTION 

Lyapunov exponents are fundamental to our understanding of chaotic processes. They 
describe the time-asymptotic rate of separation of neighboring trajectories in a dynamical 
system, and are global properties of a chaotic region, independent of the trajectory chosen to 
measure them (the fiducial trajectory). This independence is a consequence of a theorem of 
Oseledec,H which applies in the limit of infinite time. In contrast, the finite-time Lyapunov 
exponents, which describe the instantaneous rate of separation of neighboring trajectories, 
depend on the choice of fiducial trajectory and so characterize local properties of the chaotic 
region. Similarly, the characteristic direction of expansion or contraction associated with 
each exponent is a local property of the chaotic region; it converges exponentially to its time- 
asymptotic orient at ionQ'U — a much more rapid convergence rate than that of the Lyapunov 
exponents. 

The slow convergence of the Lyapunov exponents means that in practical situations, 
i.e., a computer or laboratory experiment, the system under study only rarely evolves long 
enough to "feel" the value of the true, infinite-time Lyapunov exponents. Rather, the rates of 
stretching and contracting of infinitesimal elements of phase-space are dominated by local 
effects, described by the finite-time Lyapunov exponents. In fluid dynamics and plasma 
physics, these rates are tied to the strong enhancement of transport observed in chaotic 
mixing.Q Thus, the spatial and temporal behavior of the finite-time Lyapunov exponents is 
worthy of study in its own right, in addition to their time-asymptotic value. 

In the present paper we expand upon the work of Tang and BoozeioH'Q on the geometrical 
properties of finite-time Lyapunov exponents. We consider an arbitrary smooth dynamical 
system, which we take to be a set of ordinary differential equations (ODEs), but could 
also be a discrete map. For such a system, the solution of the equations of motion is 
a smooth coordinate transformation (diffeomorphism) from the space of initial conditions 
(the Lagrangian coordinates, which label a trajectory) to the state of the system at a later 
time (the Eulerian coordinates). The coordinate transformation is smooth even when the 
system is chaotic, though it then becomes extremely complicated. 

To define Lyapunov exponents, it is necessary to define a norm on the space of our 
dynamical system, that is, a prescription for measuring distances. This norm is used to 
quantify the exponential separation of trajectories. The usual choice is the Euclidean norm, 
a mathematical embodiment of what we intuitively think of as distance. The Euclidean 



norm is also the simplest realization of a flat space. In essence, flatness is an intrinsic 
property of a space that reflects the commutativity of partial derivatives, a fact that is not 
so trivial in general coordinate systems. Flatness is taken for granted in most numerical 
and laboratory experiments, as reflected by the use of the Euclidean norm. When a space 
is not flat, it is said to have curvature; an example is a two-dimensional sphere embedded 
in three-dimensional space, used in some geophysical models where the atmosphere or the 
ocean is regarded as very thin. But presently we deal only with flat systems. 

A remarkable result of differential geometry is the existence and uniqueness of the Rie- 
mann curvature tensor, which vanishes if and only if the the underlying space is flat (See 
Sec. pni ). The Riemann tensor plays a central role in general relativity, where the curvature 
of space determines the gravitational forces on matter and light rays. Here, we shall only 
use curvature as a geometrical property; we make no use of general relativity. 

When the Riemann tensor vanishes in one coordinate system (usually, in Cartesian coor- 
dinates with the Euclidean norm), it also vanishes in any other coordinate system reached by 
a smooth transformation.! In particular, in a smooth dynamical system the transformation 
from Lagrangian to Eulerian coordinates leaves a vanishing Riemann tensor invariant. When 
the system is chaotic, so that trajectories exhibit exponential separation, some terms in the 
Riemann tensor exhibit exponential growth, while others decrease exponentially. Because 
the Riemann tensor vanishes identically, these terms have to balance each other; this requires 
the finite-time Lyapunov exponents and their associated characteristic directions to satisfy 
certain differential relations, or constraints, as was demonstrated for two-dimensional sys- 
tems by Tang and Boozer .El Under more general conditions and using the orthonormal basis 
form of the metric tensor ,13 we rederive their constraint and analyze its rate of convergence 
carefully. We also extend the calculation to three dimensions, confirming a constraint conjec- 
tured by Tang and Boozer ,H as well as deriving new ones. We support our theoretical claims 
with numerical simulations on two- (oscillating convection rollsE^I) and three-dimensional 
flows {ABC flowllil and the Lorenz modelt3). We also comment on the possible physical ap- 
plications of the constraints, specifically to the advection-diffusion equation and the dynamo 
problem. 

The outline of the paper is as follows: in Sec. |I| we review some basic concepts of 
differential geometry, as applied to dynamical systems. In Sec. |Tl| we focus our attention on 
the Riemann curvature tensor and its properties. Then, in Sec. |IV| we examine the nature of 
the vanishing curvature requirement in two dimensions. We obtain a constraint relating the 
derivatives of the smallest Lyapunov exponent and the stable direction vector. In Sec. |V| we 
carry out the same procedure in three dimensions, deriving three new constraints. We find 
that a combination of two of these constraints yields the same differential relation as in the 
two-dimensional case; the persistence of this constraint in three dimensions was conjectured 
by Tang and Boozer,a who were guided by numerical investigation of three-dimensional 
volume-preserving maps. Finally, in Sec. |V| we summarize our results, and discuss possible 
applications of the constraints. 

II. DIFFERENTIAL GEOMETRY AND DYNAMICAL SYSTEMS 

In this section we describe the point of view that the solution of a set of differential 
equations is a coordinate transformation from the set of initial conditions to the state of 



the system at a later time. In that sense, the time evolution is a map of the fluid domain 
onto itself. Moreover, if we assume that the dynamics are "smooth," the properties of 
the coordinate transformation are described by differential geometry. An analysis of these 
properties is the focus of this paper. 

We consider a prescribed smooth two- or three-dimensional vector field given by v(x, t). 
For example, v(x, t) could be a fluid flow obtained by a solution of the Navier-Stokes equa- 
tion, but in general the vector field does not have to correspond to such a fluid flow: it 
can represent the right-hand side of any smooth system of ordinary differential equations. 
However, because of the convenient framework provided by a fluid flow, we often refer to 
the phase space domain of x as "the fluid" and to infinitesimal elements of phase space as 
"fluid elements" . 

The Eulerian coordinates, x, denote the position of a point fixed in the "laboratory" 
frame. The trajectory of an infinitesimal element of phase space (also called a fluid element 
by analogy with the fluid-dynamical case) in Eulerian coordinates, x, satisfies 

^^^ = v(x(^,t),t), x(^,0)=^, (2.1) 

where ^ are the Lagrangian coordinates that label fluid elements: they are coordinates that 
describe a frame moving with the fluid. We have made the usual choice of taking as initial 
condition x(^, 0) = ^, which says that fluid elements are labeled by their initial position; 
thus, Eulerian and Lagrangian coordinates coincide at t = 0. 



The solution x = x(^,t) to Eq. (2J.) is the transformation from Lagrangian (^) to 
Eulerian (x) coordinates. For a chaotic flow, this transformation gets extremely complicated 
as time evolves. However, since we have assumed that v is a smooth vector field, the 
transformation is always different iable, so that the Jacobian matrix dx^ /d^^ is well-defined. 
This allows us to use the tools of differential geometry on the transformation x(^,t). 

Since our focus is on the rate of separation of neighboring trajectories, as character- 
ized by the Lyapunov exponents, we need to define a way of measuring distances in the 
Eulerian coordinates, i.e., a norm for vectors in the space. The norm we use is the 
Euclidean norm, defined by ||w||^ := ^"=iW*w*, where n is the dimension of the space. 
With this norm, the distance ds between two infinitesimally separated trajectories x(^, t) 
and x(^ + d$„ t) = x(^, t) + (ix(^, t) is 



ds^ = \\d^f = J2dx'dx' 



e 
We use the chain rule to write this distance in terms of the initial separation d^, 

r-^ f dx^ ,^.\ f dx^ 



<i''' = Y.{%'ii){%<ii')-Y.»n-iiV, 



where 



/ >. X v-^ dx^ dx^ ,„ ^. 

^'^•(^'^)-E7)7T7^- (2-2) 
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In the terminology of differential geometry, Eq. ( |2.2| ) is the fiat metric tensor 5ij of Eulerian 
space transformed to Lagrangian coordinates. The metric gij is a symmetric, positive-definite 
matrix that tells us the distance between two infinitesimally separated points in Lagrangian 
space. Because Qij is symmetric and positive-definite, it has n positive real eigenvalues, 
which we denote by A^(^,t), with corresponding orthonormal eigenvectors e^(^,t); we can 
then write the metric in the diagonal form 



n 



^.j(C,t) = 2^A^(e^).(e^)„ (2.3) 

where (e^)j is the ith component of e^. We assume without loss of generality that A^ > A^+i 
for yU = l,...,n — 1. The eigenvectors e^ define directions of initial separations for which 
neighboring fluid elements are converging or diverging, with A/ (^, t) the distance separating 
them. It is the rate of exponential growth of these distances that defines the finite-time 
Lyapunov exponents (or characteristic exponents) 

A^(C,t):=^lnA^(^,t). (2.4) 

(The factor of 1/2 enters because ds^ is the square of the distance.) In the limit as t — > cxd, 
the finite-time Lyapunov exponents converge to the true Lyapunov exponents, A^, which are 
independent of t and ^.Ii3 A flow is said to be chaotic if at least one of its Lyapunov exponents 
converges to a positive value. However, this convergence is very slow (logarithmic in time): 
in practical applications we are almost always dealing with finite-time Lyapunov exponents. 
For the remainder of this paper, "Lyapunov exponents" refers to the finite-time exponents, 
unless explicitly denoted as infinite-time. Whereas the Lyapunov exponents converge very 
slowly, in a chaotic region the characteristic directions e^(^,t) (the eigenvectors of g^j) 
converge rapidly (exponentially fast) to their time-asymptotic value, e°°(^).cl'Q 

In the chaotic case, we assume that Ai > (Ai > 1) and A„ < (A„ < 1), so that 
there is at least one exponentially expanding and one exponentially contracting direction. 
These are referred to as the unstable and stable directions, respectively. In three dimen- 
sions, the intermediate exponent A2 can have either sign. For bounded, autonomous fiows 
(where the velocity field v does not depend explicitly on time), the intermediate Lyapunov 
exponent must converge to zero in the infinite-time limit." For simplicity, we shall often 
use fj, = {u, m, s} to mean fj, = {1, 2, 3}, and write e^ = {u, m, s}. The letters {u, m, s} are 
abbreviations for "unstable" , "middle" , and "stable" directions, respectively. In two dimen- 
sions, we drop the extraneous middle direction and use /x = {u, s} to mean /i = {1, 2}, and 
write e^ = {u, §}. 

The existence of a contracting direction in a chaotic fiow follows from the assumption 
that the determinant of Qij 

n 

g := det gij = Yl A^ (2.5) 

stays bounded with time; otherwise, the volume of fiuid elements grows indefinitely, usually 
an undesirable situation from the physical standpoint. However, a bounded g is only a 



sufficient condition for the existence of a contracting direction, which often exists without 
this requirement. Note that g = 1 for incompressible flows (V • v = 0). 

In theory, the numerical procedure for finding the characteristic directions and expo- 
nents is straightforward: we integrate the ordinary differential equations (|2.1|) with initial 
conditions ^, together with the set of n^ equations on the tangent space, 

with initial conditions dx'^/dC,-' = S^j. The metric is then formed from dx'^ /dE,-' using the defi- 
nition ( |2.2|) , and its n eigenvalues A^ and eigenvectors e^ are obtained via standard numerical 
techniques. In practice, for a chaotic flow, the elements of the Jacobian matrix dx^/d^^ grow 
exponentially, and one quickly runs into numerical accuracy problems related to subtracting 
very large numbers. Many methods have been devised to circumvent this difficulty,lljO and 
we shall use the method of Greene and Kim.a'B'll3 Most methods yield a quantity that differs 
by a factor proportional to 1/t from the true finite-time Lyapunov exponent, and so need 
to be corrected.Q'Ej 

To illustrate the concept of the characteristic directions e^, let us take the two- 
dimensional, time-dependent periodic flow, 

~Tr^)Tr~rl! ijj(x,t) := Ak~^(smkx^ + ecosut coskx^)sm7cx'^. (2.7) 
ox^ ox^ J 

This flow is incompressible (V ■ v = 0, so that g = 1). It is used by many authors (see 
for example Refs. |10] and ^0|) to model a periodic set of oscillating convection rolls. The 
parameters are the amplitude of the rolls. A, the relative amplitude of the oscillations, e, 
the aspect ratio, k, and the frequency of oscillations, u. When uj = 0, the flow is steady, and 
the Lagraneian trajectories are nonchaotic. This is true in general of any two-dimensional 
steady fiow.ll3 

Figure |l] shows the field of stable directions s°° for the oscillating rolls, Eq. ( p.7| ). The 
two lines in the figure are portions of stable manifolds, obtained by integrating the set of n 
ordinary differential equations 

where r is the arc length along the manifold. In a chaotic flow, the stable manifold comes 
close to every point in an ergodic region.EJ In Fig. |^ we show a portion of stable manifold for 
the well-known ABC flow [Eq. ( |5.13| )], which we discuss further in Sec. 0. The interpretation 



of the stable manifold is as follows: if two infinitesimally close initial conditions lie on that 
manifold, their trajectories converge exponentially at a rate A^. 

Analogous to the stable manifold, there is an unstable manifold corresponding to u°° 
(and, in three dimensions, a manifold corresponding to m°° = s°° x u°°). The unstable 
manifold and its relation to mixing have been extensively described in Refs. ^ and ^; an 
important result is that material lines {i.e., a streak of dye) in a closed flow tends to trace 
out the unstable manifold. 



III. THE RIEMANN CURVATURE TENSOR 

We now apply some well-known results of differential geometry to the metric obtained 
in Sec. |I[ The aim is to write the Riemann curvature tensor, which characterizes the flatness 
of the space, in the basis where the metric tensor in Lagrangian coordinates is diagonal, as 
defined by Eq. (|2.3| ). In Sees. |^ and 0, we will use the invariance property of the Riemann 



tensor and the chaotic nature of the flow to establish constraints on Lyapunov exponents 
and characteristic directions. We only present the properties of the Riemann tensor that are 
essential to the argument. For more complete discussions and proofs, we refer the reader to 
standard textbooks, such as Refs. P,p^,|2^. 

Differential geometry tells us that if a metric describes a fiat space, then its Riemann 
curvature tensor 

n n 

Jim . pm pm |^ \ ^ pm -pi \ ^ pm p£ ("i A\ 

^ ijk ■— i- ji^k '- ki,j ^ 2-^ '- ki '- ji Z-^'- je.'- ki^ K"^-^) 

1=1 1=1 

must vanish in every coordinate system. (A subscript k following a comma denotes a deriva- 
tive with respect to the fcth coordinate.) The Christoffel symbols F involve derivatives of 
the metric, 

n 

^]k '■= \Z^9' igej,k + gikj - Qjk/) , (3.2) 



1=1 



where g is the matrix inverse of ga. The vanishing of the Riemann tensor is equivalent 
to requiring that second covariant derivatives commute in every coordinate system, a result 
that is taken for granted in fiat space because there covariant derivatives are the same as 
ordinary derivatives. Covariant derivatives can be thought of as derivatives that take into 
account a possible spatial dependence of the local coordinate frame. 

The Riemann tensor satisfies a number of symmetry properties, so that its total number 
of independent components is n^(?2^ — 1)/12. In three dimensions, it has six independent com- 
ponents, equivalent to the independent components of the Ricci tensor Rik := ^"=i Rhjk, 
which is symmetric. In two dimensions, the Riemann tensor has one independent compo- 
nent, equivalent the Ricci scalar R := Yl^k=i9^^ -^ik- In this paper we restrict ourselves to 
the two- and three-dimensional cases, so we only make use of the Ricci scalar and tensor. 

We are interested in finding the components of the Ricci tensor in a frame aligned with 
the local characteristic directions (the eigenvectors of gij, which are orthonormal) . The 
calculation of the Riemann tensor in such an orthonormal basisc3 simplifies considerably, as 
described by Wald.a We cite the important results here. 

Consider a set of n orthonormal vectors e, 



"M' 



E (^m)^ (^^)^ = ^M.' E (^^)^ (^/^). = '^^■' (3-3) 

i=l A'=l 

where (e^)* and (e^)^ are related by the metric via (e^)^ = ^ . gij (e^Y . In such a basis, the 
Ricci tensor is 



where the Ricci rotation coefficients are 

:=E(^'^)^M'^^(^'^)^- (3.5) 



^at^u 



«J 



We have that uja-fiu = —^aufi, so there are n^(n— 1)/2 independent cj's (two in two dimensions, 
nine in three). The symbol Vj denotes the covariant derivative with respect to the ith 
coordinate. For our purposes, Vj can be taken to be an ordinary derivative,E3 except when 
it appears as a divergence, where an extra factor involving the determinant of the metric is 
required: 

Ev.^- = E;^|7(v5n (3.6) 

so for incompressible flows there is no practical difference. 

The Ricci rotation coefficients carry information similar to the Christoffel symsbols F, 
but are defined in terms of the orthonormal basis vectors e^ instead of derivatives of the 
metric. 

Summing over the indices of Eq. (|3.4]), we obtain the Ricci scalar. 



^ = ^ ^MM = 2 X] Vi [(e^,)' Vj{e^y] - ^ {UJ^^^ Urrf, - ^ar,, ^rafj) , (3.7) 

where we have made use of the identity 



easily verified directly from the definition of tu, Eq. (|3.5| ). The Ricci tensor and scalar in an 
orthonormal frame, respectively Eqs. ( p.4|) and (|3^), are the quantities we use in Sees. ^ 
and to derive constraints on dynamical systems. 

Equation ( p.3| ) states that the vectors e^ are orthonormal with respect to a metric. We 



want to use the metric given by Eq. ( |2.3| ), so we take 

{e,)^ = Af (e,)., {e,r = A^^/^ {e,)\ (3.8) 

where, because the underlying space is Euclidean, we have (e^), = (e^)*. It is then clear 
that ( |3.3| ) is satisfied. The definitions ( |3.8| ) are called the covariant (subscripted) and con- 
travariant (superscripted) representations of the vectors e^. For a chaotic system, this 
definition of e^ implies that for long times the components (cs)* =: s* and {eu)i ='■ Ui are 
growing exponentially, and {es)i ='■ Si and (e^)* =: m* are decaying exponentially. It is these 
radically different time-asymptotic behaviors that allows the derivation of the constraints in 
Sees. and [V]. 

IV. CONSTRAINTS ON TWO-DIMENSIONAL SYSTEMS 



In Sec. Ill, we derived a form for the Ricci tensor expressed in a basis aligned with 
the characteristic directions e^ of the flow. As mentioned in Sec. |T|, the symmetries of 

8 



the Riemann curvature tensor imply that in two dimensions it has only one independent 
component, given by the Ricci scalar, Eq. ( p.7|) . In two dimensions, the terms in the Ricci 
scalar quadratic in uj cancel: 

/ ^ \y^ua^l ^TTfi — ^crTfi ^Tcrti) = "^221 "^221 + ^112 ^112 " ^221 1^221 " ^112 ^112 = 0. 
a,T,iJ, 

Hence, the Ricci scalar reduces to the simple form 

2 

i? = 2 J] V. [{e.yV,{e.y]. (4.1) 

As mentioned in Sec. |T|, we are assuming the underlying space is flat, so that R = always. 
Equation ( |4.1| ) is essentially the same expression as derived in Ref. |1|, but rewritten in a 
more transparent form and allowing for compressibility of the vector field v(x, t) (because 
here the derivatives V^ are covariant). Notice that the Lyapunov exponents enter the Ricci 
scalar through the definition of (co-)*, Eq. ( |3.8| ), as Ao- = exp(— Ao-t). 

As a direct demonstration of the vanishing of the Ricci scalar for the nonchaotic case, let 
us take the fiow v(xi, 0:2) = (/(X2), 0), a shear fiow of velocity f{x2) along the xi direction. 
For this special case, Eq. ( |2.1| ) can be explicitly solved, and the Lagrangian trajectories are 
given by 

The corresponding metric tensor is 



Edx^ dx^ 



dx^ dx^ _ ( 1 t/'(^. 



I 



2 



di^ d^^ \tf\i2) i + tV'(6; 



The eigenvalues and eigenvectors of g are then straightforward to calculate. Direct insertion 
into formula (|4.1| ) for the Ricci scalar in two-dimensions confirms, after a tedious calculation, 
that it does indeed vanish identically. When the fiow is not chaotic, the usefulness of Eq. ( [4.1[ ) 
is thus somewhat limited: it is simply an identity that is always satisfied. At best, verification 
that the Ricci scalar vanishes is a somewhat complicated consistency check. 

The real usefulness of Eq. ( [4.1| ) becomes more apparent when the Lagrangian trajecto- 
ries are taken to be chaotic (which in two dimensions requires a time-dependent fiow). It 
is still true that the Ricci scalar ( [4. 1] ) vanishes identically, because the coordinate transfor- 
mation x(^, t) is differentiable. However, there is information to be gained by considering 
what we call the curvature balance, that is, the details of how the vanishing of the curvature 
tensor is realized in the chaotic fiow. We now proceed to examine this balance. 

Let us write 

A.(et) = A~ + ^^^, (4.2) 

where the dimensionless function ?7^(^,t) satisfies limt_^oo''?^(^,i^)A = 0- The choice of the 
factor of 1/t is motivated bv the form for the time-evolution of the Lyapunov exponents 
derived by Goldhirsch et alB based on a direct analysis of the differential equations they 

9 



satisfy. Essentially, ?7^(^,t) is a "noise" term of relatively small amplitude: the dominant 
time-asymptotic behavior of the finite-time Lyapunov exponents is as 1/t. Note that with 
this definition VjA^ t = Vii]^. 

Having posed the form ([4.2|) for the finite-time Lyapunov exponents, we can rewrite ( |4.1|) 
as 



R = 2 J2 e" 



2A°°t 



V, 



e-'^^ie^yV.ie-^^ie^y) 



(4.3) 



Now we can clearly see the asymptotic time dependence of each term for large t: 



i? = 2e-2|^S°l*V, 



e ''"M* 



V,-(e-''"M^')l +2e+2|"^l*V. 



e-''^ s' 



V,{e-^^s^) 



(4.4) 



We have put absolute values around the exponents to exhibit the sign of the exponential, fol- 
lowing our assumption that there is an expanding (positive) exponent, A^, and a contracting 
(negative) exponent, A^. 

Since R vanishes identically, either the two terms in Eq. ( f4.4|) both vanish identically (as 
is the case for purely hyperbolic systems such as Arnold's cat mapti), or they must balance 
each other, growing at the same exponential rate 2 A: 

(4.5) 

In other words, the growth rate A is defined by the requirement that each term inside the 
braces in Eq. (|4.5|) be of order one asymptotically. Thus, for large time we must have 



V 
V,: 



e ''"M* 



V,{e-^-u^)\ ~exp(+2(|A^|+A)t), 
e"''^ s' V, (e-^^ s')] ~ exp(-2(|A^| - A) t). 



(4.6) 

(4.7) 



These growth rates apply even if y^ grows or decays exponentially in time, since because 
of Eq. ( |3.6| ) it appears both in the numerator and the denominator. If A > |A^|, then 
both (U) and (|3) grow exponentially. If A < |Af |, for t > (|Af | - A)~^ we have 



e-'''s'Vj{e-'^'s') 







(4.8) 



to exponential accuracy in t. In practice, we have found only cases with — |A2°| < A < |A^|. 
Equation ([4.8|) is thus a differential constraint equation relating s and rjs- 

Note that we do not have to wait for the Lyapunov exponent itself to have converged for 
the constraint (|4.8| ) to be satisfied: it is sufficient that the quantity |As| — A have a definite 
sign, which occurs rapidly in practice. 

Having derived Eq. ( [4.8|) , we now show that it can be reduced to a simpler constraint. 
Let 



K ■.= Vi{e-'^^s'). 



(4.9) 
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The constraint Eq. (f4.8| ) can be written 

(4.10) 



dK 




d 






-K\ 




= e-"^ s' V 


df 




df 





where f is a parameter that measures the distance along an s-hne, weighed by a density 
of exp(— r/s) [this relation is well-defined because exp(— r/^) > 0]. The ordinary differential 
equation ( |4.1CI| ) has solution 

K = , ^° , (4.11) 

1 + (f - fo)Ko ' ^ ^ 

which has a singularity at f — tq = —Kq^. Since all the quantities that make up K are 
well-defined for a given point on the s-line, we have to reject the nonvanishing solution to 



Eq. ( 4.10 ), and instead take K = 0. Hence, the constraint 

Vi(e-''^F)=0 (4.12) 

is necessary and sufficient for Eq. ( |4.8| ) to be satisfied. Equation ( [4.12| ) can also be written 

Vis'-s'ViVs = 0. (4.13) 

This is a compressible version of the constraint that was obtained in Ref. |l[ If we 
use In ^/g = A„t + Xgt, we can rewrite (|4.13|) as 

V-s + s- VA„t = 0, (4.14) 

where the V are ordinary derivatives {i.e., non-covariant) and any direct reference to com- 
pressibility effects drops out, but the constraint now takes the form of a relation between 
the stable direction, s, and the unstable stretching rate, A„. 

To illustrate the result, we use the oscillating convection rolls system, Eq. ( p.7|) . We 
consider the s-line starting at A and ending at B in Fig. |l], with the same parameter 
values. In Fig. |^ we plot the two terms of Eq. ( [4.14|) along that line for different times; it 



is clear that their sum converges rapidly to zero. Figure ^ highlights another property of 
the constraint (|4.13|) : if linit^oo VjS* exists (which it does in every case under consideration 
here), so must limt^oo -s' Vi?]^. 

To exhibit the rate of exponential convergence to zero of Eq. ([4.14|) , we consider the flow 



10 10 010 10 /\ 

X =x +sinx sinx , x = fix + u x +cosx cosx . (4-15) 

This is a two-dimensional autonomous system, so it cannot formally exhibit chaos. Never- 
theless, for adequate choice of fi and u trajectories of ( [4.15| ) grow exponentially away from 
the origin, so that Lyapunov exponents are well-defined even though they are not tied to any 
chaos in the system. One could also achieve such nonzero exponents with a linear system, 
but then the u- and s-lines would be straight, with vanishing derivatives, so that Eq. ( |4.14| ) 



is satisfied trivially. The nonlinear terms in in (|4.15|) do not contribute to V ■ v = i^. Fig- 
ure ^ shows the evolution of V • s + s ■ VA„t for parameter values /i = 1, z/ = — 1. The 
convergence rate is found to be — (A^ + \Xs\), corresponding to A = — A^ in Eq. (|4.7|). This 
implies from Eq. (|4.6| ) that Vi(exp(— //„) u*) remains of order one, as confirmed by numerical 
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calculation (Fig. ^ bottom). In truly chaotic systems, the convergence rate is not so readily 
obtained because the coefficient of the exponential has large fluctuations. 

The rate of convergence of the nonchaotic system seems to be an extreme case. The 
derivatives of the Lyapunov exponents along the unstable direction are well-behaved and 
do not grow exponentially, so from Eq. ([4.6[) we must have A = — A^. This leads to an 



extremely fast convergence rate for the constraint (|4.14|) . For chaotic systems, we flnd the 
derivatives of the Lyapunov exponents along the unstable direction grow exponentially, so 
the convergence rate of the constraint is slower (but still quite fast). 

The numerical evaluation of the Lagrangian derivatives of A^(^,t) and e^(^,t) is not 
straightforward. Finite-differencing (by starting two trajectories very close to each other) 
does not work well because of potentially very steep variations along the unstable direction. 
Also, it is difficult to guarantee that the two trajectories are in the same chaotic region 
when they are very close to region boundaries. A better method is to take the Lagrangian 
derivative of the set of equations used to flnd the Lyapunov exponents and eigenvectors 
themselves, which yields a set of ordinary differential equations for those derivatives. This 
is related to the method used by Tang and Boozer,El who used the Lagrangian derivative of 
the linearized system ( |2.6| ) directly. The method is unstable, but since the quantities under 
consideration converge quickly (the Lyapunov exponents are not sought, only their deriva- 
tives), reasonable accuracy can be achieved. In this paper, we have used a differentiated 
version of the continuous Gram-Schmidt orthonormalization method of Goldhirsch et alB 
with a stabilizing factor .ll§E3 The method of Goldhirsch et al. was derived independently by 
Greene and Kim,Q including the stabilizing factor. 

V. CONSTRAINTS ON THREE-DIMENSIONAL SYSTEMS 

In three spatial dimensions, the six independent components of the Riemann curvature 
tensor are embodied in the Ricci tensor (see Sec. |ITlD . As opposed to the two-dimensional 



case, the terms quadratic in the tu's in the Ricci tensor, Eq. (|3.4|) , do not vanish. 
The three diagonal elements of the Ricci tensor can be written as 

and the three off-diagonal elements as 

where {h^v^t) are ordered cyclically (i.e., {h^u^t) = {(-u, m, s), (m, s,^), (s,^, m)}) and the 
characteristic helicities are deflned as 

\/9 
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The main result of this section is as follows: if all the terms in Eq. ( |5.1| ) are bounded from 
above by exp(2|A^|t), then the constraints 

(5.3) 
(5.4) 
(5.5) 

hold to exponential accuracy for large times. Taking the difference between ( [5. 4) ) and 
and using vector identities yields 



u ■ V X rii — s 


VXmt = 


= 


m • V X u + s 


■VXut-- 


= 


u ■ 


V X u = 


= 



^V-(V^s) 



§■ VA,t 



0, 



(5.6) 



the same constraint as for the two-dimensional case, Eq. ( ^4.13| ). The constraint ( p.6[ ) was 
observed numerically for three-dimensional volume-preserving maps by Tang and Boozer. □ 
We thus see that Eq. ([5. 6]) is a consequence of two new, separate constraints. The third 



constraint, Eq. (|5.5| ), is different in nature than the previous ones, since it involves no 
Lyapunov exponents. 

If we "lift" a two-dimensional system to three dimensions by simply adding a third 
dimension corresponding to m = s x u = constant, the constraint ( ^.3|) is satisfied trivially 
because Am = and m is independent of position. The constraint ( p.4| ) can be equated 
to ( [4.13| ), the two-dimensional constraint, and (|5.5| ) is trivially satisfied because V x ii is in 
the m direction. Thus, as required for consistency, we recover the two-dimensional constraint 
from the three-dimensional ones. 

We now proceed with the derivation of the constraints ( |5.3| )-(0). The argument is 
similar to the two-dimensional case, but is complicated by the possibility that the terms 
in Eq. ( [5.1| ) and ( |5.2| ) grow at different rates, since there are more than two terms in the 
expression. There are also six independent components to the Ricci tensor, as opposed to 
just one in the two-dimensional case. 

Define 



n 



fiv 



-^fljl 



e'''V,{e,)k 



"^eiMUTiery^jVu 



(5.7) 



where rji, is defined in 



so that H 



iifj. 



{^ii/y/9)'Hf,f, and Hf,u 



Ar^ Ti/xi., ^ i^ V. 



With this definition we have separated the long-time behavior of l-L^y that is due to the 
infinite-time Lyapunov exponents. 

We focus first on the diagonal components, given by Eq. (|5.1|) . Let 2At-,- be the asymp- 
totic growth rate of the fastest-growing term (or terms) in ( p.ll ). (In two dimensions, there 
was only one common growth rate A for the two terms in i?, since they had to cancel.) 
Factoring out exp(2AT-T-t) from Eq. (|5.1|) , we obtain 



R^^ = (?^^^H^. 



+ 1 



■g-2(A.+A..)t(g^)*(^^^ _^^^) +e-2("'^+^-)*(e,)*7^^. 



_ p-2(A^+A^T)t 



{e^yn. 



+ 2e-2(^-+^-)* Huu n 



flU ' '-l/fl 



g-(A^+Ar-A,,+ATr)t -^ _ g-(Ap+A^-A„+Arr)t -^ 



g-2(A,,+A„-A^+Arr)i -^2 



(5.8) 
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where, by the definition of A,-,-, all the terms inside the braces are asymptotically either 
of order one or decreasing exponentially in time. We are interested in the terms that have 
a positive growth rate in the exponential; for such terms to remain of order one or less 
the coefficient of the exponential has to go to zero exponentially, as was the case in two 
dimensions. 

In each Rrr there is a term with coefficient exp(— 2(A^ + AT-,-)t): if A,-,- < |A^| then 
that coefficient grows exponentially at a rate exp(2(|A^| — A,-,-)t). This implies that the 
coefficient of the exponential decreases at least as fast as exp(— 2(|A^| — AT-^)t). Thus, from 
each diagonal components i?,-^, we get 



V,: 



-2f?s f.i 



s'iHr, 



H^ 



e-^^'-s' n^u\ < exp(-2(|Ari - A„Jt), from i?„, ; (5.9) 

e-'^^s^T^^™] <exp(-2(|Ar|-A^„)t), from i?^^ ; (5.10) 
+ 2?^^„?^„„<exp(-2(|A^|-A,,)t), from i?,, ; (5.11) 



The two constraints ( |5.9| ) and ( p.lO| ) can be directly inserted into ( |5.11| ), which then says 
that the product Hmu 'Hum decreases exponentially to zero. Since 



rLm.ti. f^ii 



K, 



where K is defined in ( [4.9| ), we can use the argument, proved in the two-dimensional case, 
that if Vi[s^ K] = 0, then K = 0. Thus, TCmu = 'Hum, and their product goes to zero, so 
both Timn and Hum rnust go to zero. From the definition ( |5.7|) , we see that we have just 
demonstrated the constraints ( |5.3D and ( p. 4] ), with the assumption that A^,- < |A^|. 

Also, in each Rrr there is a term with coefficient exp(2(A2° — AJ^ — A^ — A^T-)t). 
Since A^ > AJ^, the condition A < |A^| is sufficient for the growth rate to be positive. 
This leads directly to 

Hun < exp(-(Ar - A:^ + |Ari - A,.)t), 



which proves the last constraint, Eq. 

There could be other constraints, depending on the relative magnitude of the Lyapunov 
exponents and the sign of AJ^. However, the constraints derived here have the broadest 
applicability and are expected to be true in general. The requirement that A < |A^| has 
been found to hold by a wide margin for every dynamical system that we have tested. Given 
that the Lyapunov exponents enter the Ricci tensor through A~^, so that the largest growth 
rate is 2|A^|, it is reasonable to assume that the terms in the Ricci scalar could grow no 
faster than the timescale 2 A = 2|A^|. 

To exhibit the convergence of these quantities, we use the model system 



X 



x^ + sinx^ sinx^ sinx^. 



X 



x^ + cosx^ cosx^ sinx^. 



X 



Ax^ + /ix^ + z/x^. 
(5.12) 



These are an extension of the two-dimensional nonchaotic equations ( [4. 15] ) presented earlier. 
In the three-dimensional case, the system is generically chaotic, but can have unbounded 
trajectories for certain parameter values. If we pick an initial condition near the origin, then 
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initially the nonlinear terms are important, but become less so as the trajectory gets further 
away. Thus, Eqs. ( 5.12| ) have nontrivial characteristic manifolds near the origin but have 



a time-asymptotic form that is well-behaved (linear exponential growth). Figures ^ and |^ 
show the convergence of (|5.3|) and (p.4|) for different parameter values, one set incompressible 
(z/ = 0) and the other compressible (i/ 7^ 0); Figures |^ and || show the convergence of ( p.5| ) 
for the same parameter values. For this simple system, the rate of convergence for all the 
constraints is found to be — 2(Au — Am)- 

A more practical example is the ABC flow, 

v(x) = A (0, sina;^, cosa;^) + B (cosx^, 0, sinx^) + C (sina;^, cosx^, 0) (5.13) 

a sum of three Beltrami waves which satisfy V x v oc v. This flow is time-independent, 
incompressible {g = 1), and is well-studied in the context of dynamo theory; we use the 
parameter values A = 5, B = C = 2, which are known to exhibit a vigorous dynamo .liil Ell 
(All the results presented here also apply to the more common A = B = C = 1 flow, though 
the convergence rate of the constraints is slower because that flow is less chaotic.) Figure |^ 
shows the convergence of the two terms of ( ^.31) and ( p.4|) ; they are seen two track each 
other as time evolves. Figure ^ is a logarithmic plot of Hum and TCmu as a function of 
time, showing that the two quantities have converged to order 10~^. The convergence is not 
linear: it exhibits large fluctuations, yet the decreasing trend is evident. 

Finally, no demonstration of a method would be complete without applying it to the 
celebrated Lorenz model,t3 

x'^ = a(x'^ — x'^), x^ = rx^ — x^ — x^ x^^ x^ = — 6x^ + x"^ x^, (5-14) 

which is a compressible vector fleld with V-v = —(a + 1 + 6). We use Lorenz's parameter 
values of cr = 10, 6 = 8/3, r = 28. The initial condition we used is ^ = (0,1,0), which 
has a trajectory that initially skirts an unstable flxed point and undergoes nearly-periodic 
oscillations around it. Thus, the convergence of the constraints is momentarily delayed until 
the trajectory can explore the attractor, which is seen to happen at t ~ 20 in Figures ^^ 
Til. Note that this does not mean that if the trajectory comes near a flxed point again the 



convergence to zero will be undone: it will pause momentarily, to resume after the trajectory 
leaves the vicinity of the flxed point. 

The off-diagonal terms ( |5.2| ) involve the same quantities as the diagonal ones, so they 
do not immediately yield any new constraints. However, they may be useful in trying to 
establish a detailed curvature balance for the Ricci tensor (see below). 

It should be noted in closing this section that for the three-dimensional case we have not 
made use of the fact that the curvature actually vanishes. We have simply assumed that its 
individual terms grow no faster than exp(2|As|t). In contrast to the two-dimensional case, 
we have made no detailed assessment of the balance of curvature, a difficult task because 
of the number of terms involved. However, hints of this balance can often be seen in the 
numerical results, such as whenjTi^m and Timu track each other precisely in Figs. ^ and |T^, 



or in Fig. 11, where Timm and Tigs also track each other exactly. 
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VI. DISCUSSION 

Using the requirement that, in a flat space, the Riemannian curvature must vanish in any 
coordinate system, we have derived the constraints (|4.13|) in two dimensions, a generahzation 
of the work of Tang and Boozeilil to compressible flows, with a more detailed examination of 
the underlying assumptions. We have also verified the constraint for flows, whereas previous 
work had focused on area-preserving maps. Equation ([4.13|) is a constraint relating the 
stable direction, s, and the finite-time Lyapunov exponent, A^; it is satisfied to exponential 
accuracy in time in a chaotic flow. 

With the same assumptions as in two dimensions, we have derived the same constraint 
in three dimensions, Eq. (|5.6| ), verifying the conjecture of Tang and BoozeiH inspired by nu- 
merical results on volume-preserving maps. In three dimensions, it is a consequence of two 



separate constraints, Eqs. ( |5.3| ) and ( |5.4| ). There is also a third constraint, Eq. ( |5.5| ), involv- 
ing only the unstable direction, u, and no finite-time Lyapunov exponents. All have been 
shown numerically to hold with exponential accuracy in time for ABC flow and the Lorenz 
model, and we have verified the constraint for several other flows and maps, parameter val- 
ues, and initial conditions. Two of the three-dimensional constraints are satisfied trivially 
in two dimensions, and the third is just the two-dimensional constraint ( [4.13|) , showing the 
consistency of the results. 



The constraint i f).Q\) actually applies in any number of dimensions, as long as the usual 



assumption holds: all the terms in the Ricci scalar ( ^.71 ) must be bounded from above 
by exp(2|As|t). There is then a V^s^VjS-'] term in the Ricci scalar that yields the con- 
straint ( 15.61) . The higher-dimensional generalization of Eqs. (|5.3|) , (|5.4|) , and ( ^.5|) involve 
Lie brackets instead of curls, and a detailed analysis of their form remains to be done. 

Though the constraints are interesting from the mathematical standpoint, it is natural 
to ask the following: what physical irnpact do the constraints have? Previous work has 
focussed on the constraint (|5.6|) in twoa'Q and threeH'Ba dimensions. We review this work 
and discuss the consequences of the new, three-dimensional constraints. 

In Ref. |l|, it was found that the advection-diffusion equation in Lagrangian coordinates is 
a diffusion equation with an anisotropic diffusion tensor D^^ = D g^\ where D is the scalar 
molecular diffusivity and g^^ is the inverse of the metric tensor given by Eq. ( p.2|) . The dom- 
inant direction of diffusion is along the s line, with SiD^^ Sj = Dexp[—2Xst), which grows 
exponentially because As < 0. In contrast, the diffusion along the u line is exponentially 
damped at a rate exp(— 2A„t). Thus, after a short time essentially all the diffusion occurs 
along the s line. This counterintuitive result is a consequence of the diffusion equation 
smoothing out strong gradients, which are being created along the s line because of the con- 
traction of fluid elements along that direction. Thus, it is the stable direction and A^ — and 
not A^j — that characterize the chaotic enhancement to diffusion. This is where the constraint 
comes in: in Ref. |T] is was shown that, in two dimensions, Eq. ( [4. 131 ) can be written 



s ■ V{\st - In ||av||) + s ■ (k X V X K)/||Kf = 



where k = {s ■ V)s is the curvature of the § line (unrelated to the Riemannian curvature of 
the space). This implies that, along an § line, the variations in A^ are related to the variations 
in ||k;||. Indeed, this is confirmed by calculations using maps in twou and threecl dimensions: 
regions of the § line with large curvature have locally small Lyapunov exponents. We have 
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performed a similar calculation on the oscillating convection rolls system to demonstrate 
this correlation: Figure |15| (top) shows the magnitude of the curvature k as a function of 
the distance along the s line labeled A to -B in Fig. |I|, parametrized by r. The bottom part of 



Fig. |T5|is a plot of the instantaneous Lyapunov exponent A = \Xs\ (the flow is incompressible, 
so the two exponents differ only by a sign) as a function or r. Notice the local minima in A 
whenever the curvature is maximal. These locally small values of the Lyapunov exponent 
hinder the enhancement of diffusion in those regions. Such considerations were used to 
provide design criteria for a chemical reactor in Ref. ^ Numerically, the correlation was also 
verified for three-dimensional volume-preserving maps in Ref. ^ 

The constraint ( [4.13| ) has also been used in the study of the dynamo problem,c3 where a 



magnetic field is embedded in a conducting chaotic flow. Boozert^ and Boozer and Tanj 
have derived criteria that determine the spatial distribution of the magnetic field in terms 
of the magnitudes of the Lyapunov exponents. 

The constraint on the unstable direction u, Eq. (p.5|), potentially modifies these results. 



In Refs. 28 and 29, an asymptotic expression for the magnitude of the induced current in 



Lagrangian coordinates was derived, 

f ~ A3 {(B ■ u) u • V X u}' + 0(A2 A„), (6.1) 

where B is the magnetic field, and only the fastest-growing term, proportional to A^, was 
kept. The assumption was that there was nothing special about u ■ V x u, so it should be 
roughly of order one. However, we have found that typically u ■ V x u ~ A~^ A^ —>■ (see 
Figures |^, |, |Tl], and 0). Thus, the first term in Eq. ( p.lD is actually of lower order than the 



neglected term: for a nonideal plasma with resistivity rj, the power rjj'^ needed to sustain 
the dynamo is potentially much smaller than previously thought. The constraint ( |5.5| ) also 
modifies the growth rate of other quantities, such as the parallel current j ■ B. A more 
detailed analysis remains to be done, where the constraints ( |5.3D and ( p.4| ) might be found 
useful. 
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FIGURES 

FIG. 1. The field of characteristic directions s°° for oscillating rolls, Eq. ( |2.7] ), with 
A = k = e = uj = l. Two typical portions of the stable manifold are shown. 

FIG. 2. Typical portion of a stable manifold, or s°°-line, for the ABC flow with A = 5, 
B = C = 2. The domain is 27r-periodic in all three directions, as delimited by the dashed lines. 



FIG. 3. For the system of oscillating convection rolls, Eq. ( |2.7D , the convergence of the con- 
straint V • s + (s • V)Aut ^ is illustrated by evaluating its two terms on an s°°-line. Here r is 
the distance along the s°°-line from ^ to -B in Fig. IT], the solid line is — V • s°°, and the dashed 
line is (s°° • V)Aut. Some parts of the line converge more slowly because the finite-time Lyapunov 



exponent is smaller in those regions (see Fig. 15). 



FIG. 4. Top graph: Illustration of the exponential convergence of V-s + s-VA^i, for 



the nonchaotic two-dimensional model given by Eq. ( 4.15 ). The parameter values /U = 1 
and V = —1 are chosen such that the trajectories are exponentially diverging (the ini- 
tial conditions are ^ = (0.01,0.015)). The dashed line denotes an exponential with decay 
rate -(An + \\s\) ^ -(1.618 + 0.618) = -2.236. Bottom graph: plot of V • u + u • VA^ t, show- 
ing that asymptotically it remains of order one. V • u actually converges to a nonzero constant; 
the oscillations are due to u • VA, t. 



FIG. 5. TCum (solid line) and TCmu (dotted line) for Eqs. ( 5.12| ), with parameter val 



ues A = —0.1833, // = 0.75, z^ = and initial conditions ^ = (0.01, 0.015, —0.018). For comparison, 
the dashed line decreases as exp(— 2(An — Xm)t) — exp(— 2(0.698 — 0.271)t) = exp(— 0.854t). 



FIG. 6. Tlum and TCmu (respectively the solid and dotted lines, which lie on top of 
each other) for Eqs. ( 5.12|) , with parameter values A = 0.1833, // = 0.75, u = —2.6667 



and initial conditions ^ = (0.01,0.015,-0.018). For comparison, the dashed line decreases 
as exp(-2(An - A^)*) ^ exp(-2(0.396 + 0.159)t) = exp(-l.llt). The noisiness after t ~ 22 is 
due to roundoff error. 

FIG. 7. TCuu (solid line) for Eqs. ( |5.12| ), with the same parameter values and initial conditions 
as in Fig. |5|. The dashed line has the same slope as in Fig. ^. 

FIG. 8. TCuu (solid line) for Eqs. ( |5.12| ), with the same parameter values and initial conditions 
as in Fig. |6|. The dashed line has the same slope as in Fig. ^. 

FIG. 9. Top: Plot of u • V x m (solid line) and s • VA^i (dotted line) versus time for A = 5, 
B = C = 2 flow [Eqs. ( ^.13 )], with initial conditions ^ = (0.5,0.3,0.3232). The convergence of the 



two lines verifies the constraint ( |5.3D . Bottom: Plot of m • V x u (solid line) and — s • VA^i (dotted 
line) versus time for the same flow. The convergence of the two lines verifies the constraint (^.4|). 
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FIG. 10. Log-linear (base 10) plot of TCum (solid line) and TCmu (dotted line) versus time 
for A = 5, B = C = 2 flow [Eqs. ( |5.13| )], with initial conditions as in Fig. ^. There are large 
fluctuations, but the trend is for both terms to go to zero; they have reached values of 10~^ 
at t = 12. Note that, from Fig. ^, the terms that cancel to give exponential convergence to zero 
are roughly of order one, so they agree to about seven digits. For comparison, the dashed line 
decreases like exp(— 2(A„ — \m)t) — exp(— 2t). 



FIG. 11. Top: Plot of Tiuu (solid line), TCmm (dashed line), Tiss (dotted line, lying on top 
of Ti.mm) versus time for A = 5, B = C = 2 flow [Eqs. ( p.l3| )], with initial conditions as in 
Fig. |9|. For this particular system, with these parameters, we have TCmm = Tigs, though this does 
not appear to be generic. Bottom: Log-linear (base 10) plot of TCuu = u • V x u, showing its 
exponential convergence to zero — as predicted by the constraint ( ^.51 ). For comparison, the dashed 
line decreases like exp(— 2(Au — Xm)t) ^ exp(— 2t). 

FIG. 12. Top: Plot of u • V X m (solid line) and s • VA^t (dotted line) versus time for 
the Lorenz model [Eqs. ( |5.14| )], with parameter values a = 10, b = 8/3, r = 28, and initial 
conditions ^ = (0, 1,0). Bottom: Plot of m • V x u (solid line) and — s • VA^t (dotted line) versus 
time for the same flow. In both plots the two lines initially seem to converge but then diverge 
abruptly before converging again at t ~ 21. This is due to the trajectory approaching an unstable 
fixed point of the model and undergoing nearly-periodic motion. During that time the trajectory 
does not "feel" the chaos, and only after it leaves the vicinity of the fixed point does convergence 
of Tium amd Timu resume, (see Fig. p^.) 

FIG. 13. Log-linear (base 10) plot of TCum (solid line) and TCmu (dotted line) versus time for 



the Lorenz model with the same parameter values as in Fig. 12. As mentioned in Fig. 12, the 
trajectory is initially nearly-periodic. It has to first get away from the unstable fixed point and 
explore the attractor before the convergence really sets in, at t ~ 20. 

FIG. 14. Top: Plot of TCuu (solid line) and TCmm (dashed line) versus time for the Lorenz model 
with the same parameter values as in Fig. |l2|. For this particular system, with these parameters, 
we have TCmm = TCss, though this does not appear to be generic. Bottom: Log-linear (base 
10) plot of TCuu = u • V X u, showing its exponential convergence to zero — as predicted by the 



constraint (^.51). See the caption to Figs. 12 and |l^ and the text for a discussion on the transient 
oscillations before convergence sets in. 

FIG. 15. Top: For the system of oscillating convection rolls, a plot of the magnitude of the 
curvature of the s line, k = (s • V)s, as a function of the distance r along the s line labeled A to B 
in Fig. H (with the same parameter values). Bottom: Plot of the finite-time Lyapunov exponent A 
at t = 7 for the same system, as a function of r. Note the correlation between the locally small 
values of A and the magnitude of the curvature k. 
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